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Abstract 

We construct a quintom dark energy model with two non-minimally coupled scalar 
fields, one quintessence and the other phantom field, confined on the warped DGP brane. 
We study some important issues such as phantom divide line crossing, existence of the 
bouncing solutions and the stability of the solutions in this framework. We show that this 
model accounts for crossing of the phantom divide line and realization of the bouncing 
solutions. This model allows for stability of the solutions in separate regions of the uj-lo' 
phase-plane. 
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1 Introduction 



Despite all of its successes, the standard model of cosmology suffers from a series of problems. 
The most serious of these problems is the problem of initial singularity because the laws of 
physics break down at the singularity point. In order to avoid this lawlessness, there is a huge 
interest in the solutions that do not display divergencies. These solutions could be obtained 
at a classical level or by quantum modifications. Most of the efforts in quantum gravity is 
devoted to reveal the nature of the initial singularity and to understand the origin of matter, 
non-gravitational fields, and the very nature of the spacctimc. In recent analysis done within 
the loop quantum cosmology, the Big Bang singularity is replaced by a quantum Big Bounce 
with finite energy density of matter. This scenario has strong quantum effects at the Planck 
scale too. Another motivation to remove the initial singularity is the initial value problem. A 
sound gravitational theory needs to have a well posed Cauchy problem. Due to the fact that 
the gravitational field diverges at the singularity, we could not have a well formulated Cauchy 
problem as we cannot set the initial values at that point. 

On the other hand, one of the most important discoveries over the past few years is that we 
live in a positively accelerated universe which is almost spatially flat [1]. Another remarkable 
hint of the cosmological observations is that the equation of state parameter (cu) transits from 
a; > —1 to a; < —1 [2,3,4]. These discoveries generated renewed interest in bouncing models of 
the universe because it can be shown that at a positively accelerated universe a necessary condi- 
tion for a bounce in general relativity is to violate the null energy condition, i.e to have p+p < 0. 
To interpret the cosmic acceleration, a so-called dark energy component has been proposed. 
On the other hand the nature of dark energy is ambiguous. The simplest candidate of dark 
energy is a cosmological constant with the equation of state parameter to — —1. However, this 
scenario suffers from serious problems like a huge flne tuning and the coincidence problem [5]. 
Alternative models of dark energy suggest a dynamical form of dark energy, which is often real- 
ized by one or two scalar fields. In this respect, dark energy components such as quintessence, 
k-essence, chaplygin gas, phantom and quintom fields have been studied extensively [6] ( see 
also [7] and [8]). Another alternative approach to explain the universe's late-time acceleration 
is modifying the General Relativity itself [9]. Also, some braneworld scenarios are other suc- 
cessful models to achieve this goal [10]. In a braneworld scenario, our 3-brane is embedded in a 
higher dimensional bulk. Matter fields are confined to a four dimensional brane while gravity 
and possibly non-standard matter fields are free to propagate in entire space time. Among the 
braneworld models, the Randall-Sundrum 11 (RSll) model is very popular since it has a new 
modification of the gravitational potential in the very early stages of the universe evolution 
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[11]. On the other hand, the Dvah-Gabadadze and Porrati (DGP) braneworld scenario is a 
very interesting model which can describe the origin of the late-time accelerating behavior of 
the universe without adopting any additional mechanism [12]. In this setup, gravity is modified 
at large distances because of an induced four-dimensional Ricci scalar on the brane. This term 
can be obtained by the quantum interaction between the matter confined on the brane and the 
bulk gravitons. The DGP braneworld scenario explains accelerated expansion of the universe 
via leakage of gravity to extra dimension without need to introduce a dark energy component. 
While the RSII model produces ultra-violet modification to the General Relativity, the DGP 
model leads to infra-red modification. By considering the effect of an induced gravity term 
as a quantum correction in RSII model, we have a combined model that dubbed warped DGP 
braneworld in the literature [13]. This setup gives also a self-accelerating phase in the brane 
cosmological evolution. 

While DGP-inspired models essentially have the capability to explain late-time acceleration, 
crossing of the cosmological constant line and issues such as realization of bouncing solutions 
and their stability need additional mechanism to be explained in these models. With this 
viewpoint, in this paper we construct a quintom dark energy model with two scalar fields non- 
minimally coupled to induced gravity on the warped DGP brane. Wc study some currently 
important cosmological issues such as phantom divide line crossing, avoiding singularities by 
realization of the bouncing solutions and stability of these solutions. Wc analyze the parameter 
space of the model numerically and wc show that this model allows for stability of the solutions 
in the separate regions of the u)-u>' phase-plane. 



2 Warped DGP Brane 



The action of the warped DGP model can be written as follows 



S — Shuik + «5) 



^brane ; 



5 = / d'X^^g 

J hulk 



2kI 



+ 



brane 



4 



+ ^brane{gal3, "0) 



(1) 

(2) 



Here Sbuik is the action of the bulk, Sbrane is the action of the brane and S is the total action. 

with A — 0, 1, 2, 3, 5 are coordinates in the bulk while with (x = 0,1,2,3 are induced 
coordinates on the brane. k| is the 5-dimensional gravitational constant. ^^^R and ^^^Cm are the 
5-dimensional Ricci scalar and the matter Lagrangian respectively. K"^ is trace of the extrinsic 
curvature on either side of the brane. Cbraneigap, V') is the effective 4-dimensional Lagrangian on 
the brane. The action S is actually a combination of the Randall- Sundrum II and DGP model. 
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In other words, an induced curvature term is appeared on the brane in the Randall- Sundrum 
II model, hence the name Warped DGP Braneworld [13]. Now consider the brane Lagrangian 
as follows 

,,2 

(3) 



A* 



where is a mass parameter, R is the Ricci scalar of the brane, A is the tension of the brane and 
is the Lagrangian of the other matters localized on the brane. We assume that bulk contains 
only a cosmological constant, A. With these choices, action (1) gives either a generahzed DGP 
or a generahzed RS II model: it gives DGP model if A = and ^^^A = 0, and gives RS II model 
a fj, — [13]. The generalized Priedmann equation on the brane is as follows [13] 



3/1^ 



p + Po(l + eA{p,i 



(4) 



where e — ±1 is corresponding to two possible branches of the solutions ( two possible embed- 

1/2 



ding of the brane in the AdSs bulk) in this warped DGP model and A — 



where = 



l-2r) 



1/2 



P- 



.,,2£o\ 



Tj = with < 77 < 1 and po = m\ + 6^. By definition. 



POfJ- 
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mx — A^/^ and — k^"^^^- £0 is an integration constant and corresponding term in the gen- 
eralized Priedmann equation is called dark radiation term. We neglect dark radiation term in 
forthcoming arguments. In this case, generalized Priedmann equation (4) attains the following 
form, 

2?7/9\l/2 



^2 I _ 1 

?ip? 



Po 



(5) 



where p is the total energy density, including energy densities of the scalar fields and dust 
matter on the brane: 

P = + Pa + Pdm- (6) 

In what follows, we construct a quintom dark energy model on the warped DGP brane. 



3 A Quintom Dark Energy Model on the Warped DGP 
Brane 

As a part of matter fields localized on the brane, we consider a quitom field non-minimally 
coupled to induced gravity on the warped DGP brane. The action of this non-minimally 
coupled quintom field is given by 

Squint = [ d'xV^l - l^Riif'' + (7^) - Id^^di^if + Id^ad^a - Vi^p, a)] , (7) 

J brane Z Z Z 
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where ^ is a non-minimal coupling and R is induced Ricci scalar on the brane. </? is a normal ( 
canonical or quintessence) component while cr is a phantom field. We have assumed a conformal 
coupling of the scalar fields and induced gravity. Variation of the action with respect to each 
scalar field gives the equation of motion of that scalar field 

dV 

(p + + iR^ + — = 0, (8) 

and 

dv 

■& + ?>Ha - CR(J - ^ = 0. (9) 

da ^ ^ 

The energy density and pressure of the quintom field are given by the following relations 

Pquint = P,, + Pa = ^(0' - + V{ip, a) + 6CH{ipip + ud) + 3ei/'(<^' + (T^) (10) 



and 



Pquint ^ Pip + Pa ^ ^{'P'^ - o"^) " V{(f, a) - 2^{(fi(f + 2(fiH(p + (f^ + a(T + 2aH& + o"^) 



-^{2H + 3H^){^^ + a^) (11) 

In what follows, by comparing the modified Priedmann equation in the warped DGP braneworld 
with the standard Priedmann equation, we deduce the equation of state of the dark energy 
component. This is reasonable since all observed features of dark energy are essentially derivable 
in general relativity ( see [14] and references therein). The standard Priedmann equation in 
four dimensions is written as 

H^ + -^^^(pdm + Pde), (12) 

where pdm is the dust matter density, while pde is dark energy density. Comparing this equation 
with equation (5), we deduce 

Pde^ P<p + Pa + Po + spo(Al + 2ri—)\ (13) 

^ Po^ 

The conservation of the quintom field effective energy density can be stated as 

^^^g^ + 3H{pg^int + Pquint) = (14) 

Since the dust matter obeys the continuity equation and the Bianchi identity keeps valid, total 
energy density satisfies the continuity equation. In order to solve the Priedmann equation (5) 
we choose the following potential 

V{^,a)^{Cwr + lrn'{^'-<^% (15) 
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where C is a dimensionless constant describing the interaction between the scalar fields. With 
this potential, a possible solution of our basic equations, (5), (9) and (10) with supplemented 
equations (11) and (12) is as follows ( see [15] for a similar argument) 

ip — \[Cq cos(mt), a = ^JCq sin(mt) (16) 

where Co is a parameter with the dimension of mass squared describing the oscillating amplitude 
of the fields. For a fiat spatial geometry on the brane and setting pdm = 0, if we consider low- 
energy limit where by assumption p^e <^ po, we find 



(- 



Using (17) in (11), we find 



""f ~i {p^ + p„){l + ^) + po{l + eAo) 



(17) 



P^ + Pa = ^^sin\2mt) + S^H^Cq. (18) 



Therefore, Friedmann equation (17) can be rewritten as follows 

/^sin'(2mt)(l + ^) + , , , xi/2 
H^±(^ ^ 3Mf(i±£Mj . (19) 

There are four possible combinations of signs in this equation. We use this result in our 
forthcoming arguments. Before proceeding further, we note that one could choose the quantities 
in the square root in such a way that lead to a imaginary Hubble parameter. We avoid such 
cases in what follows. Also singularity points of H are treated in forthcoming arguments. 



3.1 Bouncing behavior of the model 

We start with a detailed examination of the necessary conditions required for a successful 
bounce. During the contracting phase, the scale factor a(t) is decreasing, i.e. a{t) < 0, and 
in the expanding phase we have a{t) > 0. At the bouncing point, d{t) = 0, and around this 
point d{t) > for a period of time [15,16]. Equivalently in the bouncing cosmology, the Hubble 
parameter H runs across zero from H < to H > and -ff = at the bouncing point. A 
successful bounce requires that around this point the following relation should be satisfied 

H ^ -4:nGp{l+uj) >Q (20) 

So, at the bouncing point the scale factor reaches a non-zero minimum value while the Hubble 
parameter reaches zero. By solving the Friedmann equation (19) we plot the behavior of the 
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scale factor versus the cosmic time, t, for two branches of the solutions. Figure (la) shows 
the behavior of a{t) for e = +1 and figure (16) shows the case for e = —1. As one can see, 
in both branches of this DGP-inspired model, the scale factor reaches a non-zero minimum 
and the universe switches between expanding and contracting phases alternatively. As we have 
emphasized, equation (19) has four alternative representations corresponding to four possible 
combinations of the signs. If we integrate this equation, we find 



a{t) = ttQ exp 



± 



0sin2(2mt)(l + ^) + 



PO 



3At2(l+eAo) 



1/2- 



:2i] 



Other possible combinations of signs lead to only a shift in the corresponding figures. 




0.4-1 • 1 • 1 • 1 • 1 • 1 -I • 1 • 1 • 1 • 1 • 1 • 1 • 1 • 1 

1 2 3 4 5 012345678 



(a) (b) 



Figure 1: The evolution of the scale factor for two branches of the warped DGP model with quintom 
field localized on the brane: a) Self-accelerating branch of the model ( with e = -|-1), the universe 
undergoes an expansion, reaches to a maximum radius and then crunches to a finite minimum size and 
this cycle repeats. There is no bounce at the minimum point since the scale factor has no derivative 
at that point, b) Normal branch of the model ( with e = —1). The universe switches alternatively 
between expanding and contracting phases. The minimum points are the bouncing points. 
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3.2 Crossing of the phantom divide line 

In the DGP scenario if we use a single scalar field (ordinary or phantom ) on the brane, we 
can show that the equation of state parameter of dark energy crosses the phantom divide line 
[17] ( see also [14] and [18]). It has been shown that DGP model with a quintom dark energy 
fluid in the bulk or brane, accounts for the phantom divide line crossing too [19]. Now we try 
to realize this crossing in the warped DGP braneworld with quintom matter localized on the 
brane and non-minimally coupled to induced gravity. In this warped DGP model, the equation 
of state parameter, uj of dark energy component has the following form ( with pdm = 0) 



a; = -IH ^ 



X 



Pd, 



;e 



After substituting corresponding relations for (p, a, H and V in this equation, we plot the 
behavior of oj for two branches of the DGP-inspired model versus the cosmic time. Figure 2 
shows the variation of u versus cosmic time for two possible branches of the model. In figure 
2a which is devoted to self-accelerating branch, the equation of state parameter crosses the 
cosmological constant line. This behavior is repeated periodically due to oscillating nature of 
the cosmic expansion. Figure 26 shows the situation for normal ( non-self accelerating) branch. 
In this case crossing of the cosmological constant line occurs too, but the behavior of equation 
of state parameter differs considerably compared to self- accelerating branch. As this figure 
shows, at the bouncing point uj approaches the negative infinity. 
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Figure 2: Time evolution of the equation of state parameter uj. There is a crossing of the cosmological 
constant line in both branches of the scenario: a) Self-accelerating branch, uj mimics the oscillating 
nature of the cosmic expansion, b) Normal branch. As usual, at the bouncing point uj approaches the 
negative infinity. 
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3.3 Stability of the model 



Now we study the stability of our model. The sound speed expresses the phase velocity of the 
inhomogeneous perturbations of the quintom field. In order to study the classical stability of 
our model, we analyze the behavior of the model in the u — u' plane where u' is the derivative 
of LO with respect to the logarithm of the scale factor ( see [20-23] for a similar analysis for other 
interesting cases) 

(23) 



We define the function Ca as 



(ilna dt d\a.a 



CO 

H' 



(24) 



If the matter is considered as a perfect fluid, this function would be the adiabatic sound speed 
of this fluid. But, for our model with two scalar fields, this is not actually a sound speed. 
Nevertheless, we demand that > in order to avoid the big rip singularity at the end of the 
universe evolution. Prom equation (14) we have 



Using equation of state Pde = i^dePde, we get 

Pde — ^dePde + ^dePde- 

So, the function could be rewritten as 

9 (^de 

^de- 



In this situation, the uj — cu' plane is divided into four regions defined as follows 



(25) 



(26) 



(27) 



/ : 
II : 
/// : 
IV : 



OJde > -1, 
<^<ie > -1, 
(^de < -1, 
i^de < -1, 



u' > 3uj{l + u) 
u' < 3a;(l + w) 
uj' > 3uj{l + uj) 
oj' < 3oj{l + uj) 



cl>0 
cl<0 
cl<0 
c'>0 



(28) 



As one can see from these relations, the regions I and IV have the classical stability in our 
model. We plot the behavior of the model in the oj — oj' phase plane and identify the regions 
mentioned above in figure 3. 
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-2 -1.8 -1.6 -1.4 -1.2 -1 -0.8 -0.6 -0.4 -0.2 

CO 

Figure 3: Bounds on uj' as a function oiuj in w — phase plane. The stable regions are I and IV. 

4 Summary 

One of the most serious shortcomings of the standard model of cosmology is the problem 
of initial ( and possibly final) singularity. In recent analysis done within the loop quantum 
cosmology, the Big Bang singularity is replaced by a quantum Big Bounce with finite energy 
density of matter. Also incorporation of the Gauss-Bonnet term in the action of braneworld 
models with induced gravity provides a phenomenologically rich framework to overcome initial 
singularity with possible realization of bouncing solutions [24]. On the other hand, a sound 
gravitational theory needs also to have a well posed Cauchy problem. A Model universe which 
reahzes bouncing solution is a good candidate to overcome these singularities. 

An alternative approach to explain late-time positively accelerated expansion of the universe 
is a multi-component dark energy with at least one non-canonical phantom field. The analysis of 
the properties of dark energy from recent observations mildly favor models where ^ = ^ crosses 
the phantom divide line, uj = —1 in the near past. In this respect, construction of theoretical 
frameworks with potential to describe this positively accelerated expansion and crossing of the 
phantom divide line by the equation of state parameter is an interesting challenge. In this 
paper, we have considered a quintom field non-minimally coupled to induced gravity on the 
warped DGP braneworld as a dark energy component. We have studied the bouncing behavior 
of the solutions in both branches of this DGP-inspired scenario. In the self-accelerating branch 
of the model ( with e = +1), the universe undergoes an expansion, reaches to a maximum 
radius and then crunches to a finite minimum size and this cycle repeats. In this case there 
is no bounce at the minimum point since the scale factor has no derivative at that point. In 
the normal ( non-self accelerating ) branch of the model ( with e = —1), the universe switches 
alternatively between expanding and contracting phases. The minimum points of the scale 
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factor versus cosmic time are the bouncing points. In fact there is a sequence of phases as: 
Expansion — > Turn-around — > Contraction — > Bounce and this cycle repeats regularly. 

Next we study the dynamics of the equation of state parameter. One can see that there is 
a crossing of the phantom divide line in both branches of this DGP-inspired model although 
the evolution of the equation of state parameter is different in these two branches. We have 
studied the stability of this model. As a result, there are appropriate regions oi cu — cu' phase 
plane that solutions are stable. 

Finally wc should stress on the ghost instabilities present in the self-accelerating branch 
of this DGP-inspired model. The self-accelerating branch of the DGP model contains a ghost 
at the linearized level [25]. Since the ghost carries negative energy density, it leads to the 
instability of the spacetime. The presence of the ghost can be attributed to the infinite volume 
of the extra-dimension in DGP setup. When there are ghosts instabilities in self-accelerating 
branch, it is natural to ask what are the results of solutions decay. As a possible answer we 
can state that since the normal branch solutions are ghost-free, one can think that the self- 
accelerating solutions may decay into the normal branch solutions. In fact for a given brane 
tension, the Hubble parameter in the self-accelerating universe is larger than that of the normal 
branch solutions. Then it is possible to have nucleation of bubbles of the normal branch in 
the environment of the self-accelerating branch solution. This is similar to the false vacuum 
decay in de Sitter space. However, there are arguments against this kind of reasoning which 
suggest that the self-accelerating branch does not decay into the normal branch by forming 
normal branch bubbles ( see [25] for more details). It was also shown that the introduction 
of Gauss- Bonnet term in the bulk does not help to overcome this problem [26]. In fact, it is 
still unclear what is the end state of the ghost instability in self- accelerated branch of DGP 
inspired setups. On the other hand, non-minimal coupling of scalar field and induced gravity 
in our setup provides a new degree of freedom which requires special fine tuning and this my 
provide a suitable basis to treat ghost instability. It seems that in our model this additional 
degree of freedom has the capability to provide the background for a more reliable solution to 
ghost instability due to wider parameter space. 
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